Gravitational radiation from crystalline color-superconducting hybrid stars 
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The interiors of high mass compact (neutron) stars may contain deconfined quark matter in a 
crystalline color superconducting (CCS) state. On a basis of microscopic nuclear and quark matter 
equations of states we explore the internal structure of such stars in general relativity. We find 
that their stable sequence harbors CCS quark cores with masses Af cor c < (0.78 — 0.82)Afo and radii 
fl^orc < 7 km. The CCS quark matter can support nonaxisymmetric deformations, because of its 
finite shear modulus, and can generate gravitational radiation at twice the rotation frequency of 
the star. Assuming that the CCS core is maximally strained we compute the maximal quadrupole 
moment it can sustain. The characteristic strain of gravitational wave emission ho predicted by 
our models are compared to the upper limits obtained by the LIGO and GEO 600 detectors. The 
upper limits are consistent with the breaking strain of CCS matter a < 10 -4 and large pairing 
gaps A ~ 50 MeV, or, alternatively, with a ~ 10 -3 and small pairing gaps A ~ 15 MeV. An 
observationally determined value of the characteristic strain ho can pin down the product a A 2 . On 
the theoretical side a better understanding of the breaking strain of CCS matter will be needed to 
predict reliably the level of the deformation of CCS quark core from first principles. 
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I. INTRODUCTION 

Compact (neutron) stars are important sources of 
gravitational wave radiation. Gravity waves, to lowest 
order, arise from time-dependent quadrupole deforma- 
tions of masses. Rotating, isolated compact stars will 
emit gravitational radiation if their mass distribution is 
nonaxisymmetric with respect to their rotation axis. The 
axial symmetry can be broken over detectable time-scales 
in a number of ways, e. through distortions in the 
solid phases of the star [H, 12, H, 0, [a, [E Bl, d eformations 
caused by strong magnetic fields [8l l9l. lid [Til H3 |. or 
by precession [lj, 0, EE EE S3- Continuous gravita- 
tional waves emitted by nonaxisymmetric rotating com- 
pact stars are expected to be in the bandwidth of current 
gravitational wave interferometric detectors. Upper lim- 
its on the strain of gravitational waves from a large selec- 
tion of known radio pulsars were set recently by the LIGO 
Collaboration under the assumption that the radiation is 
at twice the pulsar spin frequency [Hj]. Ellipticities de- 
rived from these upper limits are within a range that 
is compatible with theoretical predictions. Current data 
from the fifth LIGO science run places upper limits on 
the gravitational wave amplitude from the Crab pulsar 
at h = 3.4 x 1CT 25 (95% confidence), which is below the 
spin-down limit (i.e., the limit obtained assuming that 
all the energy loss is due to gravitational radiation) by a 
factor of 4 [lj| . 

In this work we study gravitational wave emis- 
sion from compact stars featuring a crystalline color- 
superconducting (CCS) phase of quark matter in their 
interiors. Our study is based on microscopic equations 
of state of nuclear and quark matter, described earlier 
in a study of the integral parameters of nonrotating and 
rapidly rotating hybrid configurations [2(j. The elliptic- 
ity of CCS quark stars and the associated gravitational 
wave emission was estimated earlier for uniform-density 
incompressible models. On the basis of such a model 



Lin [21J] concluded that the upper limit for the Crab pul- 
sar provides useful constraints on the QCD parameters. 
Haskell et al. [22| estimated core deformations and as- 
sociated with them constraints on the QCD parameters 
for sequences of 1.4 M Q and R = 10 km stars containing 
CCS quark cores with different transition densities from 
nuclear to quark matter. The nuclear matter equation of 
state was approximated by a n — 1 polytrope, and quark 
matter was treated as an incompressible fluid. Equilib- 
rium and stability of hybrid compact stars constructed 
from microscopic equations of state of quark and nuclear 
matter were studied in Ref. [20l |. For stellar sequences 
that are based on microscopic input the stability is not 
a guaranteed feature. The equilibrium among the quark 
and nuclear phases is achieved for rather stiff nuclear 
equations of state, while those configurations that con- 
tain quark phase(s) often belong to unstable branches of 
the sequences [II (H H [H, [27], H, [H . For our micro- 
scopic input a stable branch of CCS hybrid stars emerges 
in the form of a "second family" of configurations, which 
is separated from their purely nuclear counterparts by 
a region of instability. The masses of the resulting hy- 
brid configurations are close to the maximum sustainable 
mass M max ~ 2M Q , i.e., are much more massive than 
"canonical" neutron stars with M ~ 1.4M Q . 

Below, by resorting to microphysical input we add real- 
ism to the treatment of the gravitational radiation from 
hybrid compact stars, in particular, the constraints on 
microscopic parameters derived from the observational 
upper limits. Even though we lift a large portion of un- 
certainties related to the microphysical input of the the- 
ory, the elastic properties of the CCS phase and other 
solid phases of our models (nuclear crusts) remain uncer- 
tain. The breaking strain that solid phases of compact 
stars can sustain is difficult to compute and is usually as- 
sumed to lie in the range 1(T 5 < er < 1(T 2 . The values 
close to the upper "optimistic" bound of cr implies distor- 
tions in the star's crusts and/or core that are interesting 



from the perspective of gravitational wave emission. We 
concentrate below on distortions in the CCS phase and 
will ignore any other contributions, as those originating 
from the distortions in the crusts, precession or magnetic 
deformations. Distortions may exist in other types of 
compact objects, such as solid strange stars, quark stars, 
stars with mixed phases of quarks, meson-condensates 
and hadrons [30|, |3l|, |32|, |33[ ; the problem of gravitational 
wave emission from such objects is discussed in the lit- 
erature [13, [35| . Furthermore, we assume that the CCS 
quark matter is the true ground state of quark matter 
in some density interval [H [13, H [H, So, HO, El, El 
and that the value of the dynamical strange quark mass 
is small en oug h to favor the three-flavor variant of this 
phase [H |40L ItO. El. However, among the candidate 
phases for the ground state of quark matter at interme- 
diate densities, the three-flavor CCS quark phase and its 
two- flavor counterpart are the only phases that behave as 
solids (nonzero shear modulus). Therefore, their possible 
manifestations via gravitational wave emission is unique 
among the color-superconducting phases. 

This paper is organized as follows: In Sec. HI] we study 
the equilibrium sequences of nonrotating hybrid com- 
pact stars featuring CCS quark cores. We compute the 
masses, radii, and quadrupole moments of the quark 
cores of hybrid configurations. Section IIIII discusses the 
gravitational wave radiation from the CCS quark cores. 
The experimental upper limits on the strain of gravita- 
tional waves from pulsars are compared to the theoretical 
predictions of our models. Our conclusions are collected 
in Sec. HVl 



II. MODELS OF CCS QUARK STARS 

We consider models of hybrid compact stars con- 
structed from zero-temperature microscopic equations of 
state of nuclear matter described by Dirac-Bruckner- 
Hartree-Fock (DBHF) theory and quark matter de- 
scribed by the Nambu-Jona-Lasinio model (see also 
Ref. [III). The DBHF approach is based on a self- 
consistent solution of the equation for the in-medium rel- 
ativistic T matrix and nucleon self-energy starting from a 
realistic nucleon-nucleon potential. The potential is ad- 
justed to reproduce the nuclear phase shifts and deuteron 
binding energy. The saturation properties of nuclear 
matter are well reproduced within DBHF theories. The 
quark equation of state is derived from the Lagrangian 
of the Nambu-Jona-Lasinio model. This model is based 
on the four-quark contact interaction picture, i.e., the 
gluons are integrated out from the theory. It exhibits 
the chiral symmetry restoration in quark matter, but 
lacks confinement. Its coupling constants and momen- 
tum cutoff are adjusted to reproduce the meson masses 
and the pion decay constant in the vacuum. Dcconfine- 
ment phase transition between the nuclear and quark 
phases is implemented via a Maxwell construction, i.e., 
it occurs at the point where the pressures of both phases 
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FIG. 1: (Color online) Number density versus pressure for 
models A [solid, black online), Al [dashed, red online), and B 
(dashed-dotted, blue online). For models A and Al the nuclear 
(low density) equation of state is the same; for models A and 
B the quark (high-density) equation of state is the same. At 
the deconfinement phase transition there is a jump in the 
density at constant pressure. 



as functions of the baryonic chemical potential coincide. 
At the deconfinement phase transition the density expe- 
riences a jump at constant pressure and zero tempera- 
ture. Within the quark matter phase a Cooper pair con- 
densate of quarks with mismatched Fermi-surfaces will 
emerge. Our models feature the so-called three-flavor 
Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) or crystalline 
color-superconducting (CCS) phase as a possible candi- 
date of the ground state of Cooper-paired quark mat- 
ter [H, M, H El, E3, El, EI El. In this phase the 
condensate order parameter is spatially modulated; as a 
result the superfluid phase behaves as a solid with finite 
shear modulus. This, in turn, implies that the LOFF 
phase can sustain equilibrium nonaxisymmetrical defor- 
mations, which potentially can lead to gravitational wave 
radiation. It should be noted that a more complete treat- 
ment of high-density matter will need to include phases 
other than the CCS phase and the possibility of tran- 
sitions between them 



44]. However, the latter phases 



are not crystalline and are unlikely to affect the physics 
of gravitational wave emission in a direct manner. Indi- 
rectly, they may have an effect, e. g., by changing the 
fractional volume occupied by the CCS phase. 

We consider three different equations of state that are 
displayed in Fig. [T] (sec also Fig. 1 of Ref. 20]). For 
models A and Al the nuclear (low density) equation of 
state is the same; for models A and B the quark (high- 
density) equation of state is the same. At the deconfinc- 
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FIG. 2: (Color online) Dependence of the total stellar mass 
and the mass of the quark core in units of solar mass Mq 
on the central density for nonrotating configurations based 
on equations of state A, Al, and B (the labeling is as in 
Fig. [y. The lower set of curves represents the masses of 
the CCS quark cores, the upper set - the total masses of the 
configurations. The maximal masses are marked with boxes. 



merit phase transition there is a jump in the density at 
constant pressure. The high-density regime contains two 
equations of states for crystalline color superconductiv- 
ity which differ by the normalization of pressure at zero 
density (bag constant). Model Al is normalized such 
that the pressure vanishes at zero density. For models 
A and B the zero density pressure is shifted by amount 
Sp = 10 MeV/fm 3 . The quark condensate is assumed to 
be in the strong coupling regime, where the ratio of the 
couplings in the quark-quark and quark-antiquark chan- 
nels rj = 1 (the "canonical" value is rj — 0.75 and follows 
from the Fierz transformation from the quark-anti-quark 
to quark-quark channel.) The values of the (baryonic) 
chemical potential (in MeV) and pressure (in MeV/fm 3 ) 
at the quark-nuclear matter interface are (1234.2; 96) for 
model Al J1230.0; 95) for model A and (1234.8; 108) for 
model B [2(|. We start by considering the equilibrium 
and stability of nonrotating cold hybrid stars with crys- 
talline color-superconducting cores. Our emphasis will 
be on the internal structure and the integral parameters 
(mass, radius, etc) of the quark cores, rather than the 
integral parameters of the stars (these are discussed in 
Ref. [20(). We parameterize the sequences of equilibrium, 
nonrotating stellar configurations in general relativity 
in terms of the central density p c of the configuration. 
It is assumed that the configurations are cold, i.e., the 
stellar structure does not depend on temperature. (Be- 
cause of large magnitudes of gaps, color-superconducting 
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FIG. 3: (Color online) Dependence of radii of nonrotating 
compact stars and the radii of their quark cores on central 
density for models A, Al, and B (the labeling is the same as 
in Fig. [TJ . The lower set of curves represents the radii of the 
CCS quark cores, the upper set the radii of the configurations. 



phases can affect the early thermal evolution of neu- 
tron stars; we assume that our models have cooled down 
to temperatures well below the respective Fermi ener- 
gies [45|, HH, H3])- We obtain the spherically symmetric 
solutions of Einstein's equations for self-gravitating fluids 
by solving the well-known Tolman-Oppcnhcimcr-Volkoff 
equations 48]. A sequence of configurations is stable if 
the derivative dM/dp c is positive, i.e., the mass is an 
increasing function of the central density. At the point 
of instability the fundamental (pulsation) modes become 
unstable. If the stability is regained at higher central 
densities the modes by which the stars become unstable 
toward the collapse belong to higher order harmonics. 

The masses of the stellar configurations (M) and the 
CCS quark cores (M corc ) as a function of the central 
density are shown in Fig. [5] The sequences of hybrid 
configurations (A/ corc > 0) appear with increasing cen- 
tral density when the density of deconfinement phase 
transition is reached. Prior to that point the stars are 
purely nuclear and become unstable (dM/dp c < 0) be- 
fore the stability is regained by the hybrid configurations. 
The masses of this "second family" vary in a narrow 
range around the maxim attainable mass 2M Q , which is 
reached at the central density 2 x 10 15 g cm~ 3 . The val- 
ues of the maximum masses are large as a consequence 
of the hardness of the underlying nuclear equations of 
state. The masses of the CCS quark cores cover the 
range < M COTe /M Q < 0.75 — 0.88 for central densi- 
ties 1.3 x 10 15 < p c < 2 x 10 15 . Thus, the quark core 
mass ranges from one third to about the half of the total 
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FIG. 4: Dependence of density on the internal radius for 
model A. The dashed lines correspond to purely nuclear stars, 
the solid lines to hybrid configurations. For latter models 
there is a density jump at the phase transition to quark mat- 
ter. The transition radius saturates at the value R COIC = 7 
km as the central density is increased. 



stellar mass. Above one solar mass quark cores can be 
harbored only by unstable configurations. 

Fig. [3] displays the radii of hybrid configurations along 
with the radii of the CCS quark cores. With the onset of 
CCS quark phase the radii of the configurations become 
smaller, i.e., the hybrid configurations are more compact 
(and more massive) than their neighboring purely nu- 
clear counterparts. Within their stability range the hy- 
brid configurations have almost constant radii R = 12 
km. The radius of the quark core increases rapidly with 
the onset of the quark phase and saturates at the asymp- 
totic value i?corc = 7 km. The quark cores occupy about 
20% of the volume of the star. The density profiles of 
configurations are displayed in Fig. 2] For purely nuclear 
configurations the density profile is flat in the core. As 
the central density exceeds the deconfinement transition 
density, the profiles show a density jump at some internal 
radius. The jump reflects the behavior of the equation of 
state at the deconfinement phase transition. There is a 
rapid growth of the radius (volume) of the quark phase 
as the density is increased from 1.3 x 10 15 to 1.6 x 10 15 
g cm -3 . For larger central densities, i.e., more massive 
objects, the radius at which the transition from quark to 
nuclear matter takes place is independent of the central 
density of the configuration (see also Fig. [3]). 

Figure [5] displays the mass-radius {M-R) relationship 
for the stellar configurations (right set of curves) and 
the CCS quark core (the left set of curves) for models 
A, Al, and B. The M-R diagram for the stellar con- 



FIG. 5: (Color online) Mass-radius relationship for hybrid 
stars for models A, Al, and B (the labeling is as in Fig. [TJ. 
The right set corresponds to the stellar configurations, the 
left set to the CCS quark cores. 



figurations has been discussed elsewhere [2Cj; we recall 
that the obtained mass-radius relationships are consis- 
tent with the current astronomical bounds on neutron 
star masses. They are consistent with the line in the M- 
R diagram inferred from EXO0748-676 [49[. Note that 
the emission of redshifted lines used to put constraints 
in Ref. (49l | has not been confirmed in a subsequent ob- 
servation of Ref. [5(| . Also note that alternative models 
of quark stars are consistent with these bounds, which 
therefore do not rule out the hypothesis of quark mat- 
ter in compact stars [5l|. The M-R tracks for the CCS 
quark core display the rapid increase in the radius of the 
core in the presence of the small amount of quark matter 
[the M(R) slope is almost horizontal to the j?-axis]. The 
maximal masses of quark cores of the stable configura- 
tions are in the range M < 0.75 — 0.88M Q depending on 
the equation of state; cores with larger masses correspond 
to configurations that are unstable toward collapse into 
a black hole. 



III. QUADRUPOLE MOMENTS AND 
GRAVITATIONAL WAVE EMISSION 

The characteristic strain amplitude of gravitational 
waves emitted by a triaxial star rotating about its prin- 
cipal axis is 



where v is the star's rotation frequency, r is the distance 
to the observer, e = (I xx — I yy )/I zz is the equatorial 
ellipticity, /y is the tensor of moment of inertia, G is 
the gravitational constant, c is the speed of light. The 
characteristic amplitude, Eq. |T]), does not involve the 
orientation of the source with respect to the observer; it 
is related to the sources' angle averaged field strength 
(h 2 ) = Jdfl(hl + /i|)/4tt by the relation h ~ 1.15</i 2 ). 
Alternatively, the strain amplitude can be expressed in 
terms of the m = 2 mass quadrupole moment as 



h 



16tt 2 G /32tt 
15" 



1/2 



Q 22 ^ 2 



(2) 



where 
J 



the 



mass multipole is defined as Qi m = 
J pr l Y l l n d 3 r. Thus, to compute the characteristic strain 
amplitude of gravitational waves we need to compute the 
quadrupole moment of a CCS quark core in static equilib- 
rium. The elastic deformations are assumed to be small 
perturbation on the background equilibrium of the star; 
the equilibrium between gravity and elastic forces, in the 
Cowling approximation, implies that 0] 



V^Ty = 5pg(r)fj, 



(3) 



where g(r) is the local gravitational acceleration, 5 rep- 
resents Eulerian perturbation, = —pgij + Uj is the 
stress-energy tensor of the CCS phase, tjj is the shear 
stress tensor, and <7y is the flat 3-metric (indices i and 
j run over 1,2,3). Upon expanding the tensor t t j in 
spherical harmonics and integrating ([3]) with the bound- 
ary condition that assumes vanishing stress in the nuclear 
envelope surrounding the CCS quark core one finds [cf. 
Ref. 0, Eq. (64)] 
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where [7 = 2 + d\ng(r) / dhir and t rri t\ and t r ± are the 
coefficients of the expansion of the shear stress tensor 
in spherical harmonics 0]. These (unknown) coefficients 
can be determined if one assumes that the CCS core is 
maximally strained. The shear modulus is defined as 
(i = tij/2<7ij, where <7ij is the strain tensor. The shear 
modulus of CCS matter has been estimated in Ref. [42j 

M = 2.47 MeV far 3 ( A ) ' f ^ ) 2 (5) 

' \10 MeV J V400 MeV/ ' 1 ; 

where A is the gap parameter, p q is the quark chemical 
potential. For a maximally strained CCS core 
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FIG. 6: (Color online) The quadrupole moments of stellar 
configurations computed according to Eqs. (T3J), © and © 
for models A, Al, and B as a function of central density (the 
labeling of models is as in Fig. [!}. Each triple of curves cor- 
responds to the gap parameter value (from top to bottom) 
A = 50, 40, 30, 20, and 10 MeV. The magnitude of the break- 
ing strain is <r max = 10~ 2 , and Q ma x scales linearly with er ma x- 



where <7 ma x is the maximal value of the quantity a 2 = 
(TijO*3 /2 (breaking strain). It is assumed that the strain 
in the CCS core is position independent. By combining 
Eq. ^ and ^ we obtain the maximal quadrupole mo- 
ment Q max . If the core is uniform and incompressible, 
then 13511 



Qr 



13^0-max-Rco 
G-Mcore 



(7) 



Figure [6] displays maximal quadrupole moments of a 
sequences of hybrid stars derived by combining Eqs. |4|. 
([5]) and The quadrupole moment depends on the 

product of the breaking strain and the shear modulus. 
The former will be treated as a free parameter in view 
of the poor knowledge of this quantity; in Fig. [S] we as- 
sume <7 max = 1CP 2 . Since the equation of state fixes the 
quark chemical potential, the only parameter on which 
the shear modulus depends is the gap, which we will 
vary in the range 10 < A < 50 MeV. It is assumed that 
the gap is independent of the density of the CCS quark 
phase. Treatments of the CCS phase that are based on 
the Ginzburg-Landau expansion near the critical tem- 
perature predict values of gaps that are in the range 
5 < A < 25 MeV. The low temperature CCS phase (more 
relevant for the problem at hand) is likely to support a 
larger, up to 100 MeV, gap. The maximal quadrupole 
moment scales as Q max ~ CT max A 2 according to Eqs. (@|)- 
©■ 




FIG. 7: The quadrupole moments of stellar configurations 
(solid lines) computed for fixed values of the shear modu- 
lus of CCS quark matter fi = 5 and fi — 50 MeV fm~ 3 for 
model A as a function of central density. The dashed lines 
represent the uniform-density incompressible matter approx- 
imation, Eq. (0, where the masses and radii are set equal to 
that of the quark cores of microscopic models (the values of 
the central density on the x axis do not represent the central 
densities of these models) . 



Figure [Jj compares the quadrupole moments of quark 
cores of stellar configurations computed from the mi- 
croscopic equations of state A, Al, and B with the 
quadrupole moments of uniform, incompressible objects 
computed via Eq. ((7|) . The masses and radii of the latter 
are set equal to that of the quark cores in the microscopic 
models. (Note that the central densities on the abscissa 
of Fig. [Jj refer only to the microscopic models.) It is seen 
that the uniform, incompressible approximation system- 
atically overestimates the magnitudes of the quadrupole 
moments, although numerically the error is not large. An 
accurate value of the quark core radius, that follows from 
our microscopic treatment, appears to be a more impor- 
tant factor in providing realistic values of Q m ax because 

Of Qmax ~ #c re Scaling. 

Figure [5] displays the strain of gravitational wave emis- 
sion given by Eq. ([2|) for models A, Al, and B as a func- 
tion of the core mass M corc (the results for model A over- 
lap with that of model B on the scale of the figure) . Each 
pair of curves corresponds to the gap parameter value 
(from top to bottom) A = 50, 40, 30, 20, and 10 MeV. 
We assume rotation at the frequency of the Crab pul- 
sar v — 29.6 Hz at a distance of 2 kpc. Since ho is a 
linear function of quadrupole moment Q maX j its depen- 
dence on the gap function reflects the quadratic depen- 
dence of Qmax on the gap function discussed above. The 



FIG. 8: (Color online) The strain of gravitational wave 
emission for a pulsar rotating at the Crab pulsar frequency 
v — 29.6 Hz at a distance of 2 kpc, computed according to 
Eqs. (UJ) and ([5]) for models A and Al (the results for models 
Al and B overlap on the scale of this figure). Each triple of 
curves corresponds to the gap parameter value (from top to 
bottom) A = 50, 40, 30, 20, and 10 MeV. The magnitude 
of the breaking strain is <r max = 10 -4 , and ho scales linearly 
with o" max . The current upper limit on the strain of the Crab 
pulsar derived by the LIGO experiment is shown by the hor- 
izontal dashed line. 

magnitude of the breaking strain is <7 max = 10 -4 . The 
current upper limit on the characteristic strain of gravi- 
tational radiation from the Crab pulsar is 3.4 x 10~ 25 and 
is shown in Fig.[8]by the horizontal dashed line. It is seen 
that the strain of gravitational wave emission is close to 
the upper limit for A = 50 MeV. If, however, the gaps are 
small more "optimistic" values of 10 ~ 3 < a max < 10~ 2 
will be needed to generate sizeable strain amplitude. 
Clearly, because of the large uncertainty in the actual 
value of (T m ax no definite conclusion can be drawn about 
the microscopic parameters of the CCS phase from the 
current upper limits on the strain of gravitational wave 
emission. The current upper limit for gravitational waves 
from the Crab pulsar implies cr max A 2 ~ 0.25 MeV 2 (un- 
der the assumptions of the present model) . 

Figure [9] shows the strain of gravitational wave emis- 
sion ho for two approximations to the quadrupole de- 
formations (realistic vs incompressible), discussed above 
and illustrated in Fig. [Jj As in Fig. [7j the masses and 
radii of incompressible models are set equal to that of the 
quark cores derived from the microscopic models, i.e., the 
abscissa of Fig. [9] refers to central density of the micro- 
scopic models only. It is seen that the uniform, incom- 
pressible approximation systematically overestimates the 
magnitude of ho for all central densities. 
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FIG. 9: Same as in Fig. [8] but for fixed values of the shear 
modulus of CCS quark matter /i = 5 and pi = 50 MeV fm -3 
for models A (solid lines) and Al (dashed lines). The re- 
sults for models Al and B overlap on the scale of the figure. 
For each pair of curves with the same value of /j, and same 
equation of state the lower curve is obtained from the full 
expression for Q ma x, while the upper one from the uniform- 
density incompressible fluid approximation, Eq. (J7J. The 
latter approximation systematically overestimates the value 
of the quadrupole moment. The magnitude of the breaking 
strain is a" max = 10 -4 ; ho scales linearly with a ma x- 



IV. CONCLUSIONS 

We studied gravitational wave emission by hybrid 
compact (neutron) stars harboring maximally strained 
crystalline color-superconducting quark cores on a ba- 
sis of microscopic equations of state. In a first step, 
we (re)constructed hybrid configurations, which are com- 
posed of color superconducting quark matter at high and 
purely nuclear matter at low densities. The high-density 
quark matter is described in terms of the semimicroscopic 
Nambu-Jona-Lasinio model which includes pair correla- 
tions that lead to the three-flavor LOFF phase as the 
ground state of QCD at moderate densities. The low 
density nuclear phase is described in terms of hard rela- 
tivistic equations of state based on the Dirac-Bruckner- 
Hartree-Fock theory. These configurations in static grav- 



itational equilibrium arise as a "second family" of stable 
configurations. We find that the stable (with respect to 
a collapse into a black hole) segment of the sequences 
have CCS quark cores with masses < M covc /Mq < 
0.75-0.88 for central densities 1.3 x 10 15 < p c <2x 10 15 
and radii up to 7 km. Consequently, about 1/3 of the 
mass of hybrid configuration is concentrated in its quark 
core, which occupies about 20 % of star's volume. 

In the second step, we considered the characteristic 
strain of gravitational wave emission due to deformations 
of the CCS quark core. If the quark core is maximally 
strained, the components of the stress tensor can be ex- 
pressed through product of the breaking strain and the 
shear modulus of the CCS quark matter. The quadrupole 
moment and the characteristic strain mainly depend on 
the not-well-known product a max A 2 . The current upper 
limits imply that <r max A 2 < 0.25 MeV 2 , which can be ac- 
commodated by assuming either large gaps A ~ 50 MeV 
and moderate values of <7 max ~ 10~ 4 or small gaps ~ 15 
MeV and more "optimistic" values of breaking strain, 
ff max ~ 10 . In view of large uncertainty in the value of 
breaking strain 10~ 5 < cx max < 10~ 2 no definite conclu- 
sions can be made about the possible magnitude of the 
gap in the CCS phase. However, if the gaps are large, 
A > 50 MeV, the values cx max > 10~ 3 are excluded by the 
present scenario of gravitational wave radiation. Note, 
however, that the assumption that the CCS core is max- 
imally strained may not correspond to the actual situa- 
tion in compact stars and the evolutionary avenues that 
may lead to such maximal deformations are not explored 
yet. One may only speculate that a rapidly spinning pul- 
sar at birth may be triaxially stretched, a deformation 
that may be preserved within the star after the transi- 
tion to the crystalline state. Nevertheless, we conclude 
that the characteristic strain amplitudes of gravitational 
waves produced by our models are well within the reach 
of the current gravitational wave detectors for a wide and 
reasonable range of parameters. 
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